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Abstract Roche-lobe overflow and common envelope 
evolution are very important in binary evolution, which 
is believed to be the main evolutionary channel to hot 
subdwarf stars. The details of these processes are diffi- 
cult to model, but adiabatic expansion provides an ex- 
cellent approximation to the structure of a donor star 
undergoing dynamical timescale mass transfer. We can 
use this model to study the responses of stars of various 
masses and evolutionary stages as potential donor stars, 
with the urgent goal of obtaining more accurate stabil- 
ity criteria for dynamical mass transfer in binary popu- 
lation synthesis studies. As examples, we describe here 
several models with the initial masses equal to 1 Mq 
and 10 Mq, and identify potential limitations to the 
use of our results for giant branch stars. 
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1 Introduction 

The Fourth Meeting on Hot Subdwarf Stars and Re- 
lated Objects was convened in Shanghai, China from 
19th to 24th, July. Recent discoveries and develop- 
ments in both theory and observation of hot subdwarfs 
and related objects were reported, and many unsolved 
problems were discussed. As we all know, hot subd- 
warf stars are extreme horizontal-branch stars or re- 
lated objects. They may dominate the UV-upturn of 
early-type galaxies and they exist in both the field of 
our Galaxy and its globular cluster syste m. About half 



of subdwarf B (sdB) st ars are binaries (JMaxted et al. 



2001t ISaffer et al.l 120011 ) ; binary evolution is obviously 
important in their formation, as the observed systems 
are too compact to have avoided pa st mass transfer . 



Binary populat ion synthesis models (IHan et al.l 12002 . 
2003tlHanll2008 ) explain naturally the sdB binary frac- 



tions and the UV-upturn of early-type galaxies via: (a) 
one or two phases of common envelope (CE) evolu- 
tion, (b) stable Roche-lobe overflow (RLOF), and (c) 
the merger of two He-core white dwarf stars (WDs). 

RLOF and CE evolution are also very important in 
the formation of other binary systems, e.g., cataclysmic 
variables, X-ray binaries, double white dwarfs and bi- 
nary neutron stars. Those systems containing compact 
objects are among the most energetic and rapidly vari- 
able sources known. Unfortunately, we know very little 
about the details of CE evolution, but dynamically un- 
stable RLOF appears to be the trigger that launches 
it. The threshold conditions for dynamical mass trans- 
fer now in common use in binary population synthesis 
calculations are b ased on polytropic models for rapid 
mass loss process ( Hjellming fc Webbinklll987L and ref- 
erences therein). These polytropic studies provide use- 
ful qualitative insights into RLOF and CE evolution, 
but they omit much relevant physics, and fail to ad- 
dress many advanced evolutionary stages of interest. 



Detailed studies o f binary evolution (|Han et alj 12002 : 
Chen fc Hanll2008l ) reveal a need for more realistic de- 
terminations of the threshold conditions for dynamical 
mass transfer. 

With the motivation and cautions above, we set out 
to study st ellar rapid mass loss based on the pioneer- 
ing work of iHiellmina (|1989a. bl. We describe the ba- 
sic assumptions and numerical techniques employed in 
modeling stellar adiabatic mass loss in section [5J with 
initial results and their possible application to binary 
population synthesis in section [3l In section IH we offer 
a short discussion of some remaining problems. 



2 Stellar adiabatic mass loss model 

Dynamical time scale mass transfer occurs when, in re- 
sponse to mass loss, the interior expansion of a donor 
star drives its surface beyond its Roche lobe. That ex- 
pansion is driven by the quasistatic readjustment of the 
star to the local drop in pressure as overlying mass is 
removed, and it is characterized by the star's dynam- 
ical time scale. That time scale is typically many or- 
ders of magnitude shorter than either the nuclear or 
thermal time scale of the donor star. In dynamical 
time scale mass transfer, therefore, the response of the 
donor star to mass loss becomes asymptotically adia- 
batic |Ge et al.ll20"l0l : iHiellming fc Webbinklll987t ). and 



its composition profile remains fixed. Inside its Roche 
lobe, the donor remains in hydrostatic equilibrium, ex- 
cept very near the inner Lagrangian point. Indeed, even 
for thermal timescale mass transfer, thermal relaxation 
is largely confined to the outermost layers of the donor, 
and those layers are quickly stripped away. Model se- 
quences in which potential donor stars are stripped of 
mass and allowed to respond only adiabatically there- 
fore provide valuable insight into the inner structure 
of donor stars undergoing rapid mass transfer. Fur- 
thermore, we can determine whether a given donor star 
will be stable or unstable to dynamical time scale mass 
transfer, by comparing its radius as it loses mass adia- 
batically to the effective radius of its Roche lobe. 

The stellar structure equations for adiabatic mass 
loss models are time- independent. The equations for 
hydrostatic equilibrium and mass continuity are re- 
tained, but those for energy conservation and energy 
transport are replaced by algebraic constraints fixing 
the entropy and composition profiles: 



dlnP 
dm 

dm 



Gm 

' A-KT'^P' 



2Trrp' 



(1) 



(2) 



s(rn) — so(rn), 



X(m) = Xo(m), 



(3) 



(4) 



with structure variables P for pressure, r for radius, 
m, for mass. T for temperature, s for specific entropy 
and X for the abundances of various nuclear species. 
The initial specific entropy sq, and abundances Xq are 
tabulated as functions of the radial mass coordinate to 
in the initial model for an adiabatic mass loss sequence, 
and remain unchanged. 

The solution to these equations not only describes 
the dynamical response of a star to mass loss, but it pro- 
vides valuable insights into how that star then relaxes 
thermally in response to adiabatic expansion. Because 
the specific entropy and composition are fixed as func- 
tions of mass, so too arc their gradients with respect to 
mass. One can write the local energy fiux through the 
stellar interior in terms of the local entropy gradient, 
and thereby reconstruct the interior luminosity profile. 
Gradients in that luminosity profile in turn imply a spe- 
cific disposition of energy sources and sinks through the 
stellar interior. 



3 Initial results 

Based on the stellar evolu tion code de v elope d by 
EggletonI |l97ll . Il972l . Il973[ ) and IPaxtonI (|2004h . we 
built a numerical code in FORTRAN95 to calculate 
stellar adiabatic mass loss models. The latest input 
physics of our co d e are t he same as in Eg gleton's code 



( Hanet all Il994 12003: iPols et al.l Il99 



i ll998D . We 
set the mixing-length parameter of a = l/Hp = 2.0 
and the convecti ve overshooting parameter S^v = 0.12 
( Pols et al.lll998[ ). With our stellar adiabatic mass loss 
code, we plan to study the responses of donor stars to 
mass loss at different evolutionary stages and metallic- 
ities of interest. As an initial application, we describe 
here two sets of models with metallicity Z = 0.02, ini- 
tial masses 1 Mq and 10 Mq, respectively. 

3.1 1 Mq star: radiative core with a convective 
envelope 

In its simplest terms, the adiabatic response of a star 
to mass loss depends on the star's specific entropy pro- 
file. In radiatively stable zones, entropy increases out- 
wards, and the removal of high-entropy surface layers 
brings lower-entropy (denser) material toward the sur- 
face, leading to contraction; in convective zones, the 
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entropy profile is nearly flat (except possibly for su- 
peradiabatic surface layers) and the star expands in re- 
sponse to mass loss. Stars typically consist of both 
radiative and convective layers, so there is a competi- 
tion between these effects, but one should note that, 
in strictly adiabatic mass loss, the entropy profile re- 
mains fixed in mass, and so convective boundaries also 
remain fixed in mass. With these concepts in mind, let 
us examine the behaviour of a 1 Mq star in different 
evolutionary stages as it loses mass. Such a star has 
a surface convection zone throughout its lifetime, but 
that zone contains relatively little mass on the main se- 
quence, deepening dramatically as the star reaches the 
giant branch. Fig. [T] shows the response of a 1 Mq star 
to mass loss at different evolutionary stages: 

1. Surface convection zones are invariably capped by 
a superadiabatic layer in which convection becomes 
very inefficient. Until the star begins ascending the 
giant branch, however, the mass in this superadia- 
batic zone is quite small (about 10~^ Mq); the star 
brightens as higher-entropy material is exposed, but 
its change in radius is negligible (see (a), and (b) in 
Fig.[l]). 

2. ZAMS: Chemically- homogeneous lower main se- 
quence stars are reasonabl y well-described by a 



composite polytrope model (JRappaport et al.lll983 : 



Hjellming fc Webbinklll987[ ). This 1 Mq ZAMS star 



contracts because of the sharp specific entropy gra- 
dient in the outer part of its radiative core. After 
about half of its initial mass is lost, however, the en- 
tropy profile of the remaining star has become very 
fiat (see solid line in Fig. [5]), and the star begins to 
expand in response to mass loss. 

3. ZAMS to TMS: As the star evolves through the main 
sequence, the specific entropy of the core decreases, 
while that of the outer envelope increases (see dashed 
line in Fig. [2]), and the contraction in response to 
mass loss becomes more pronounced. 

4. TMS to RGB: The mass and depth of the convective 
envelope increases dramatically as the star reaches 
the RGB. Its radial response becomes dominated by 
its convective envelope, and star expands modestly 
in response to mass loss. 

5. RGB to TGB: Though the mass fraction of the con- 
vective envelope becomes smaller, its specific en- 
tropy grows enormously as the star evolves; convec- 
tion becomes very inefficient near its surface, and 
the star develops an extended superadiabatic sur- 
face layer. That layer dominates the initial response 
of the star to mass loss, leading to a rapid initial ex- 
pansion. The response of the stellar radius to mass 
loss does not show a simple relation with mass frac- 
tion of the convective envelope as do main sequence 
stars. 



3.2 10 A/0 star: convective core with a radiative 
envelope 

Although the 10 Mq models have very different inner 
structure from the 1 Mq models, the behaviour of these 
models (see Fig. ^ can be understood with the same 
general precepts: 

1. The relatively massive 10 Mq stars lack a surface 
convection zone until late in their evolution. Spe- 
cific entropy increases dramatically very near the 
surface, as seen in FiglH As these layers are stripped 
away, the star decreases rapidly in effective tempera- 
ture and luminosity as lower-entropy material is ex- 
posed, but the amount of mass involved is still small 
(less than about 10~* A/©), and so the stellar radius 
changes very little at first. 

2. ZAMS to TMS: 10 Mq main-sequence stars have 
pronounced positive entropy gradients in their radia- 
tive envelopes, leading to stronger initial contraction 
than seen above in the 1 Mq models. However, they 
also have large convective cores, with fiat entropy 
profiles that ensure they re-expand with sufficient 
mass loss, like the 1 Mq models. 

3. TMS to RGB: As 10 Mq stars evolve across the HR 
diagram, the rapid rise in specific entropy through 
the envelope softens. Their initial contraction in 
response to mass loss becomes less severe, but the 
growing entropy difference between core and enve- 
lope (see Fig. |4|) suppresses re-expansion. 

4. RGB to TGB: Like 1 Mq stars, 10 Mq stars with 
deep convective envelopes expand in response to 
mass loss (compare Fig. [3] (d) and Fig. |4] (d)). At 
their high luminosities and large radii, they, too, de- 
velop superadiabatic surface layers, and expand vi- 
olently in response to its removal. 

3.3 Application to binary population synthesis 

When a star fills its Roche lobe {R ~ -Rl), the sta- 
bility of the binary against dynamical time scale mass 
transfer depends on whether adiabatic perturbations to 
the system are damped or unstable. That criterion for 
stability can be written as (|Webbinklll985l ) 



Cad > Cl, 



(5) 



where Cad = {dh\R/ d\nM)ad is the stellar adiabatic 
radius-mass exponent, and Cl = (91ni?L/i91nM)ad is 
the time-independent first-order response of the Roche 
lobe to mass transfer. In the simplest models of mass 
transfer, in which systemic mass and orbital angu- 
lar momentum loss rates are proportional to the mass 
transfer rate, Ql depends only on the binary mass ratio. 



Our stellar adiabatic mass loss models therefore allow 
us to determine the critical mass ratio (donor/accretor) 
above which a binary will be unstable to dynamical time 
scale mass transfer. This is an essential consideration in 
mapping out binary population synthesis models. We 
are currently building a libra ry of stellar adiabatic mass 
loss models ( Ge et al.ll201Cll) . covering the full range of 
evolutionary stages from ZAMS to the tip of the RGB 
or AGB, as appropriate. Fig. [5] shows the model grid 
of our project. 



4 Discussion and conclusions 

RLOF process and CE evolution are very impor- 
tant in the formation of hot subdwarf stars and 
other binary systems . Ear lier published studies of 
Hiellming fc WebbinkI (|l987l . and references therein) 
provided useful qualitative insights in RLOF and CE 
evolution, but left considerable room for improvement, 
especially with regard to the range of evolved phases 
that need to be addressed. These deficienci es were in 
part redressed in later work ( Hiellminall989b), but only 
a fragm ent of that work was ever published ( Hiellminel 
1989a[). The stellar adiabatic mass loss models de- 



scribed in this paper considerably extend the scope of 
Hjellming's work. They allow us not only to study the 
interior structure of donor stars undergoing dynami- 
cal timcscale mass transfer, but also to evaluate the 
stability criteria for dynamical mass transfer in binary 
population synthesis. 

An initial application of our models to 1 Mq and 
10 Mq stars has been described in this paper. The 
reader should beware that the these models have lim- 
itations to their usefulness, as they depend on being 
able to separate the process of dynamical relaxation to 
hydrostatic equilibrium from that of thermal relaxation 
to thermal equilibrium. On the main sequence, for ex- 
ample, the dynamical time scale of the Sun is roughly 
10"^" that of its thermal time scale, making the adi- 
abatic mass loss model an excellent description of the 
asymptotic behavior of a solar-type star to very rapid 
mass loss. However, at the opposite extreme in the HR 
diagram, stars near the tip of the RGB or AGB have 
thermal time scales approaching their dynamical time 
scales, and the entire mass transfer process may require 
full-blown time-dependent modeling. This convergence 
of thermal and dynamical time scales is broadly re- 
lated to the abrupt expansion seen in our models of 
1 Mq and 10 Mq stars at the tip of the giant branch 
(TGB), when their outermost superadiabatic layers are 
stripped away. In reality, Roche lobe overflow is far 



from spherically-symmetric, as treated here (by neces- 
sity) , but that superadiabatic expansion may neverthe- 
less reflect a real physical phenomenon. Convection it- 
self is not spherically symmetric, and with convective 
velocities in the superadiabatic zone approaching sound 
speed, it may be possible for rising flows near the inner 
Lagrangian point to bridge the potential to the compan- 
ion's Roche lobe while the donor still lies well within its 
own Roche limit. We are not now able to pass judge- 
ment on that possibility. 

When dynamical instability occurs, common enve- 
lope evolution almost certainly follows, as the thermal 
time scale of the accreting star is invariably longer than 
that of the donor, which itself is generally much longer 
than its dynamical time scale. This ordering of time 
scales ensures that the envelope cannot cool efficiently 
on the transfer time scale, but remains extended and 
engulfs both stellar cores. It may happen that thermal 
time scale mass loss from the donor is still rapid enough 
to form a common envelope, but this case can only arise 
if the system is stable against dynamical mass trans- 
fer, while unstable to thermal time scale mass transfer. 
Since the stellar dynamical time scales of both donor 
and accretor are shorter than their thermal time scales, 
the prospect arises of a quasistatic common envelope, 
that is, of formation of a contact binary. Such large 
numbers of such objects are known (the W UMa sys- 
tems, as examples) that they must be very long-lived, 
evolving in a very different fashion from CE evolution 
as we have used that term above. 

We believe that stellar adiabatic mass loss models 
provide the most useful approach to date toward defin- 
ing the limits of dynamical stability in interacting bina- 
ries, essential input to the construction of binary pop- 
ulation synthesis models. 
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Fig. 1 Hertzsprung-Russell (HR) diagrams and radial re- 
sponses of 1 Mq models during adiabatic mass loss. Frames 
(a) , (b) and (c) show the tracks in the HR diagram of stars 
initiating mass loss at the zero-age main sequence (ZAMS 
- solid line), in the Hertzsprung gap (HG - dash-dot line), 
and at the tip of the red giant branch (TGB - dash-dot- 
dot-dot line), respectively. The solid squares and numbers 
in each of these frames indicate the amount of mass lost 
(in Mq) at that point in the track. Frame (d) shows the 
radial responses to mass loss for each of the models shown 
in frames (a)-(c), with the addition of models at the ter- 
minal main sequence (TMS - dashed line) and at the base 
of the red giant branch (RGB - dotted line). Models with 
relatively shallow surface convection zones (ZAMS, TMS, 
and HG) initially contract in response to mass loss, while 
those with deep convective envelopes (RGB and TGB) ex- 
pand, dramatically so when convection becomes inefficient 
and the outer layers have strongly superadiabatic tempera- 
ture gradients. 
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Fig. 2 The specific entropy profiles of the initial 1 Mq 
models of the mass- loss sequences shown in Fig. [T] Lines are 
coded as in Fig.[T]according to evolutionary phase at the be- 
ginning of each mass loss sequence. As the star evolves, the 
specific entropy in its core decreases. A sharp positive en- 
tropy gradient develops just outside the core with the onset 
of hydrogen shell burning, becoming more pronounced with 
time. The steep negative entropy gradients at the surfaces 
of more evolved stars are due to superadiabatic (inefficient) 
convection. 
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Fig. 3 HR diagrams and tlie radial responses of 10 M© 
models during adiabatic mass loss. Frames (a), (b), and 
(c) show tracks in the HR diagram of stars initiating mass 
loss at ZAMS, HG, and TGB phases, respectively, as in 
Fig. [l] The solid squares and numbers in each of these 
frames indicate the amount of mass lost (in Af©) at that 
point in the track. Frame (d) shows the radial responses 
to mass loss for each of the models shown in frames (a)- 
(c), as well as for TMS and RGB models. Lines are coded 
according to evolutionary phase as in Fig. [T] 
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Fig. 4 The specific entropy profiles of the initial 10 Mq 
models of the mass-loss sequences shown in Fig. [S] Lines 
are coded as in Fig. [l] according to evolutionary phase at 
the beginning of each mass loss sequence. As seen in Fig. O 
the specific entropy in the core decreases as the central com- 
position changes. The flat entropy profiles in the core are 
due to convection, while the steep positive gradients at the 
surface are due to the radiative envelope. As seen in Fig. [2l 
sharp positive entropy gradients are seen outside hydrogen- 
burning shells, and in the case of the TGB model, we see 
also the sharp increase outside the helium-burning shell. 
The TGB model again has a very extended, superadiabatic 
envelope, reflected in the sharp negative entropy gradients 
near the surface. 
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Fig. 5 The model grid for different stars with metallicity 
Z = 0.02. The full set of models (not shown) ranges in 
mass from 0.1 M© to 100 Mq, at intervals of logM = 0.1. 
Solid circles, squares and triangles indicate different evolu- 
tion stages for high, intermediate and low mass stars, re- 
spectively. They mark initial models for adiabatic mass loss 
sequences. 
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